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ABSTRACT 

We study the nature of rapidly star-forming galaxies at 2; = 2 in cosmological hydro- 
dynamic simulations, and compare their properties to observations of sub-millimetre 
galaxies (SMGs). We identify simulated SMGs as the most rapidly star- forming sys- 
tems that match the observed number density of SMGs. In our models, SMGs are mas- 
sive galaxies sitting at the centres of large potential wells, being fed by smooth infall 
and gas- rich satellites at rates comparable to their star formation rates (SFR). They 
are not typically undergoing major mergers that significantly boost their quiescent 
SFR, but they still often show complex gas morphologies and kinematics. Our simu- 
lated SMGs have stellar masses of M* ~ 1O"-" '^M0, SFRs of - 180 - 500 Mg/yr, 
a clustering length of ^ 10/i~^Mpc, and solar metallicities. The SFRs are lower than 
those inferred from far-IR data by ~ x3, which we suggest may owe to one or more 
systematic effects in the SFR calibrations. SMGs at z = 2 live in ~ IO^^Mq halos, 
and by z they mostly end up as brightest group galaxies in ~ lO^^Mg halos. 
We predict that higher-Af* SMGs should have on average lower specific SFRs, less 
disturbed morphologies, and higher clustering. We also predict that deeper far-IR 
surveys will smoothly join SMGs onto the massive end of the SFR— M* relationship 
defined by lower-mass z 2 galaxies. Overall, our simulated rapid star-formers pro- 
vide as good a match to available SMG data as merger-based scenarios, offering an 
alternative scenario that emerges naturally from cosmological simulations. 

Key words: galaxies: formation, galaxies: evolution, galaxies: high redshift, galaxies: 
starburst, submillimetre, methods: N-body simulations. 



1 INTRODUCTION 



Sub-millimetre galaxies (SMGs; iBlain et al. I l2002h are 
among the most enigmatic objects in the Universe. SMGs 
were originally identified using the Submillimetre Common- 
User Bolometer Array (SCUBA) at 850/im and 450/im on 
the James Clerk Maxwell Telescope, with multi-wavelength 
follow-up from the radio to the X-rays. Their most remark- 
able property is their enormous bolometric luminosity, some- 
times exceeding W^^Lq, most of which is emitted as dust- 
reprocessed light in the far infrared (IR). Because the sub- 
millimetre band probes a falling (versus wavelength) portion 
of the galaxy's spectral energy distribution (SED), SMGs 
have the fortuitous property that their apparent brightness 
is relatively invariant over a large redshift range. In prin- 
ciple, this can be exploited to probe star formation and/or 



black hole growth over a large fraction of cosmic time, but 
in practice this requires ancillary data to both locate the 
object precisely and to obtain its redshift. Radio interfer- 
ometry has proven useful for this, enabling precise location 
of an optical counterpart for which a redshift can be ob- 
tained. In this way, the majority of SMGs are found to lie 
at 2: ~ 1.5 — 3 (|Aretxaga et al.l|2007h. with the most distant 
one currently known at 2 = 4.76 ( Coppin et al.|[2009l) . How- 
ever, the sensitivity of current radio interferometers dies out 
at 2; > 3, meaning that such a selection may bias the over- 
all redshift distribution to lower redshifts. Recent data using 
AzTEC at 1.1 suggests a somewhat higher mean redshift 
for these sources than t he 850- /^m selected ones by includ- 
ing 24/im counterparts (|Chapin et al. |[2009l ). In any case, it 
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is clear that SMGs probe the brightest galaxies during the 
most active cosmic epoch for galaxy assembly. 

Observational studies of SMGs are among the fastest- 
growing areas of extragalactic astronomy. Uniform and com- 
plete samples of SMGs are now being obtained such as 
the SCUBA Half-Degree Extragalactic Survey (SHADES; 
ICoppin et al.l 2006) covering 0.5 square degrees to a depth 
of 2 mjy, and th ere are other large samples such as that of 
IChapman et al.l (|2005 '). Such surveys can in principle pro- 
vide inter esting discriminants betw een models for the origin 
of SMGs (|van Kampen et al.ll2005l ). With new instruments 
such as LABOCA on the Atacama Pathfinder Experiment, 
AzTEC, and eventually SCUBA-2 on JCMT, and upcom- 
ing facilities in the near future such as the Large Millime- 
ter Telescope, catalogues of SMGs will continue to improve 
rapidly. 

Despite rapidly accumulating data, the nature of the 
SMGs remains poorly understood. Locally, the most bolo- 
metrically luminous galaxies are ultraluminous IR galax- 
ies with luminosities of Lir >~ IQ^^Lq (ULIRGs; 



ISanders fc Mirabel I Il996l '). and are seen to be ongo- 
ing major gas-rich mergers, e.g. Arp 220. By the First 
Law of Duck^, it was (and often continues to be) as- 
sumed that SMGs are analogous merger events at high 
redshift, simply scaled up to larger luminosities owing 
to the higher gas content of early galaxies. The ma- 
jor merger hypothesis is s upported by the pr eponderance 
of close neighbours (e.g. Ilvison et al. I 2007h and their 
often disturbed morpholo gies ( Menendez-Delmestre et al. I 
I2OO7I : iTacconi et~all I2OO8I '). Although recent models sug- 
gest a strong coevolutio n of stellar and black hole growth 
in m ajor mergers (e.g. Idi Matteo. Springel, fc Hernguist I 
I2OO5I ). observations indicate that star formation domi- 
nates the bolometric luminosity, and active galact ic nuclei 
(AGN) provide only a minor contributio n (Alexan der et al.l 
20051: iMenendez-Delrnestre et al. I [2OO7I : [Pope et al. I I2OO8I : 
Clements et al. I I2OO8I ). This isn't necessarily a difficulty 



for the merger scenario if the peak star formation occurs 
at an earlier phase in the merger relative to peak black 
hole g rowth, as such models suggest (e.g. Narave man et al.] 
l2009al '). But it does mean that the far-IR luminosity of 
SMGs can be plausibly translated into a star formation rate 
(SFR) without the fear of large contamination by AGN. If 
one adopts local calibrations (SFR= 4.5 x 10~*'*Lfir erg/s; 
iKennicutt I [l998a ) . the inferred SFRs are typically many 
hundreds to several thousands of solar masses per year. 

Hierarchical models of galaxy formation have tradition- 
ally had difficulty reproducing the observed numbers and 
fluxes of SMGs. Two theories have figured prominently in 
recent discussions, both exploiting the large enhancement 
in star formation rates believed to accompany gas-rich ma- 
jor mergers: (1) SMGs are large, gas-rich, merger-induced 
starbursts, caught in a special phase where their luminos- 
ity is significantly enhanced over their quiescent state (e.g. 
iNaravanan et al.. ,2009a ) and (2) SMGs are modest-sized 
merger-induced starbursts whose bolometric luminosity is 
greatly enhanced by an extremely top- heavy stellar initial 
mass function (IMF: lBaugh et al.|[2005h . 

The first scenario, large major mergers, is the canoni- 



cal one, and it has had recent success reproducing the de- 
tailed S EDs and CO properties of SMGs from merger sim u- 
lations (IChak rabarti et al.ll2008l : INaravanan et al." 2009a'bh. 
However, owing to the rarity of such large galaxies at high-z 
and the short duration of their merger-induced boost, this 
scenario may have difficulty reproducing the number den- 
sity of SMGs within a hierarchical context . To be quantita- 
tive, to produce SMG fluxes (~ 5 mJy), INaravanan et all 
(|2009ah needed to merge galaxies with stellar masses of 
~ 3 X IO'^^Mq. Such galaxies are above L « at these epochs 
and at z ~ 2 the stellar mass fu nction of iMarchesini et al.l 
(2009) and lKaiisawa et all (|2009l ) indicates that their num- 
ber density is « 5 x 10~^ Mpc~^ (comoving). In addi- 
tion, the dimensionless major merger rate for such ob- 
jects is around unity per Hubble tim e at z ^ 2 — 3 (even 
allowing for up to 3 : 1 mergers; IGuo fc White I I2OO8I : 
iHopk ms ilaDllooi), meaning that if such objects are ob- 
servable as SMGs for 1 00 Myr (as indicated in Figure 1 of 
INaravanan et al.ll2009al ). then Sit z — 2 the predicted num- 
ber density of SMGs would be ~ 2 x 10"'' Mpc"^. The 
observed number density is up to an order of magnitude 
higher than this , 1 - 2 x lO"'^ Mpc"^ (|Borvs et aL II2005I : 
ISwinbank et al . 2006: Dvc ct al. 2008). Similar merger simu- 
lations bv lChakrabarti et al.. (,2008i ) found that even higher 
masses are required to produce SMG fluxes, which would 
result in even lower abundances. The stellar mass functions 
at 2 ~ 2 are still s ubject to signiflcant systematics (e.g. 
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iMuzzin et"al~ll2009l: iKaiisawa et aLll2009l '). so the discrep- 
ancy could be less. iGenel et al. I (|2008l 7 found that there are 
enough halos with sufficient accretion rates undergoing ma- 
jor mergers to account for perhaps a quarter to a half of 
the observed SMGs, albeit with some assumptions about 
how r apid galaxies merge compared to their halo merger 
times. iDekel et al.l (|2009l ) have also proposed that rapidly 
star-forming galaxies at high redshift are powered primarily 
by smooth accretion and minor mergers, though they sug- 
gested that some SMGs themselves might still be powered 
by major mergers. Hence while large gas-rich mergers can 
certainly produce objects that look like SMGs, it is possible 
that a large fraction of SMGs are in fact not such objects. 

An alte r native scenario for SMGs was presented by 
iBaugh et all (120051 ). based on semi-analytic galaxy forma- 
tion models. Owing to the paucity of larg e major mergers 
in hierarchical models, iBaugh et al.l (|2005l ) were driven to 
argue that more common mergers of sub-L* galaxies give 
rise to SMGs. But such modest-sized mergers fell far short 
of reproducing the observed SMG fiuxes. To boost the flux, 
they then invoked an IMF that is flat above IMq during the 
merger event. This model can match both the number den- 
sity and far-IR fluxes of SMGs (although the actual SFRs are 
well below those inferred using standard conversion factors), 
at the cost of introducing a radically top-heavy IMF in such 
systems. However, subsequent data on SMGs has not sup- 
ported this model. For instance, this scenario predicts stel- 
lar masses signiflcantly lower t han that recently infer red for 
SMGs from near-IR data (e.g. ISwinbank et al.lliooil ') . Also, 
the IMF directly inferred from CO and dynamical observa- 
tions of SMGs disfavours such drama tic departures from a 
standard (e.g. IChabrier I l2003h IMF (|Tacconi et al.ll2008l ). 
ISwinbank et alT i 20081 ) discusses some possible ways to rec- 
oncile these discrepancies, but a fully consistent scenario has 
yet to be put forth. 



The Nature of Suh-mm Galaxies 3 



A t hird, less dis c ussed scenario for SMGs was pre- 
sented in iFardal et all l|200ll ). This work considered the na- 
ture of SMGs in cosmological hydrodynamic simulations, 
and showed that for reasonable values of dust temperatures 
such simulations can reproduce the observed 850/im num- 
ber counts. In these models, SMGs were galaxies forming 
stars at high rates, > 100 Mq /yr, but were massive systems 
doing so on relatively long timescales rather than in short 
me rger-induced bursts. A similar conclusion was reached 
by iFinlator et afl (|2006| ). In this work on the properties 

01 z = A galaxies from cosmological hydrodynamic simula- 
tions, it was noted that a few galaxies had extremely high 
star formation rates, including two exceeding 1000 Mq/jt. 
These galaxies were not undergoing a large burst of star for- 
mation owing to a major merger, but instead were massive 
galaxies that had been forming stars at hundreds of Mq/jt 
for some time. Such high SFRs are expected at high red- 
shifts because the dense intergalactic medium (IGM) and 
short cooling times yield large accretion rates ( | Keres et al.l 
I2OO5I : iDekel fc Birnboim II2OO6I : iKeres et al.ll2009l ). This sce^ 
nario, therefore, suggests that SMGs are not associated with 
major mergers at all, but are instead super-sized versions of 
normal star-forming galaxies. 

In this paper, we follow on the works of iFardal et al.l 
(|200ll ) and IFinlator et al~l (|2006l ') to examine in detail the 
nature of simulated galaxies with high SFRs at redshift 

2 ~ 2, and how their properties compare to the wealth 
of recently obtained multiwavelength data on SMGs. We 
employ cosmological hydrodynamic simulations including a 
heuristic implementation of galactic outflows that has been 
shown to match observations of more typical z ~ 2 galaxies 
reasonably well, described in ^ In ^ we study the stel- 
lar masses, star formation rates, gas fractions, metallicities, 
environments, and clustering of simulated SMGs, and show 
that the "super-sized star-former" scenario comes interest- 
ingly close to reproducing many of their observed properties. 
In ^we study four simulated SMGs in more detail, includ- 
ing their morphology, kinematics, and star formation and 
enrichment histories. A significant challenge to our scenario 
is that the SFRs inferred for SMGs are factors of few higher 
than in our simulated SMGs at a fixed number density. In [J5] 
we suggest that SMG star formation rates have been mod- 
estly overestimated, and place this idea in a broader context 
of galaxy evolution at those epochs. We summarise and dis- 
cuss the implications of our results in 



2 SIMULATIONS 

We employ our modified version of the N- 
body-|- hydrodynamic code Gadget -2, described more 
fully in Qppenheimer fc Dave I (|2008l ): here we review the 
basic ingredients. Gadget-2 uses a tree-particle-mesh 
algorithm to compute gravitational forces on a set of par- 
ticles, and an entropy-conservi ng formulation of Smoothed 
Particle Hydrodynamics (SPH: ISpringel iboOSi ) to simulate 
pressure forces and shocks on the gas particles. We include 
radia tive cooling from prim ordial (following | Katz et al] 
1 19961 ) and metal (based on [Sutherland fc Dopita I Il993l ') 
elements, assu ming ionis a tion equilibrium. Star forma- 
tion follows a JSchmidt I (|l959l ) Law calibrated to the 
iKennicutt I (|l998l ) relation; particles above a density 



threshold where sub-particle Jeans fragmentation can 
occur are randomly selected to spawn a star with half 
their original gas mass. The interstellar medium (ISM) 
is modelled through an analytic subgrid recipe following 
iMcKee fc Ostriker lll977[). including en ergv returned from 
supernovae (|Springel fc Hernguist hoO^ ). 

Chemical enrichment is followed in four individual 
species (C, O, Si, Fe) from three sources: Type II SNe, Type 
la SNe, and stellar mass loss from AGB stars. The former in- 
stantaneously enriches star-forming particles, whose metal- 
licity can then be carried into the IGM via outflows (de- 
scr ibed below) . Type la modelling is based on the fit to data 
bv lScannapieco fc Bildstenl(|2005h . as a prompt component 
tracing the star formation rate and a delayed component 
(with a 0.7 Gyr del ay) tracking stellar mass. Stellar mass 
loss is derived from iBruzual fc Ch ariot ' ("2003") population 
synthesis modelling assuming a Ichab ricr (2003) IMF. De- 
layed feedback adds energy and metals to the three nearest 
gas particles; Type la's add 10^^ ergs per SN, while AGB 
stars a dd no energy, only metals. We use lLimongi fc Chieffi I 
(200^ yields for Type II SNe; yield s for Type la and AGB 
stars come from various works (see lOppenheimer fc Dave I 
|2008| . for details). Solar metallicity is, in our case, defined as 
a total metal mass fraction of 0.0189. 

Kinetic outfiows are also included emanating from all 
galaxies. Gas particles eligible for star formation can be ran- 
domly selected to be in an outfiow, by which the particle's 
velocity is augmented by in a direction given by vxa, 
where v and a are the instantaneous velocity and accelera- 
tion. The ratio of probabilities to be in an outflow relative 
to that to form into stars is given by 77, the mass loading 
factor. Hydrodynamic forces on wind particles are "turned 
off" until a particle reaches one-tenth the threshold density 
for star formation, or a maximum time of 20 kpc/v^, which 
attempts to mock up chimneys where outflows escape that 
would otherwise be unresolvable in cosmological simulations 
of appreciable volume. 

The choices for u„ and 77 define the "wind model". 
Here we use scalings expected for mom entum-driven 
winds (|Murrav. Quatert. fc Thompson 1 12005| ). though note 
that such scalings can be generated b y other physical 
scenarios (|Dalla Vecchia fc Schavel 120081 '). In this model, 
Hiu = Say/ — 1, and rj = ao/a, where a is the veloc- 
ity dispersion of the host galaxy identified using an on- 
the-fly galaxy find er. We estimate a from the galaxy's stel- 
lar mBss_JollowingJ^to^ Mao, fc White! (ll998D . as described 
in lOppenheimer fc Dave I l|2008l ). /l is randomly chosen to 
lie betwee n 1.05 — 2, based on observations of local star- 
bursts by iRupke. Veilleux fc Sandersi l|2005h . This choice 
also gi ves outflow velocitie s consistent with that seen in 
z ^ 1 (IWeiner et al. I I2OO8I ) and 2 ~ 2 - 3 l|Steidel et al. I 
|2004 ) star- forming galaxies, cro is a free parameter, and is 
adjusted to broadly match the cosmic star formation his- 
tory; we choose ao = 150 km/s. This yields 77 ~unity for L* 
galaxies at 2 = 2, which is in ag reement wit h the constraints 
inferred from 2 ~ 2 galaxies bv lErbl l|2008h . 

This momentum-driven wind model has displayed 
repeated and often unique success at matching a range of 
cosmic star formation and enrichment data. This includes 
observations o f IGM enrichment as seen i n 2 ~ 2 — 4 
Civ systems JOppenheimer fc Dave I l2006l. I2OO8I ). 2 ~ 
O VI systems (jOppenheimer fc Dave 1 I2OO9I ) and 2 ~ 6 
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meta l-line absorbers (|Oppenheimer. Dave, fc Finlatoj 
|2009| ). Concurrently, it also matches the galaxy mass- 
meta llici ty relation |Dave. Finlator. fc Qppenheimer I 



I2OO7I : iFinlator fc Dave 
richment levels 



200^ . 



and 



the en- 

intragroup 

gas (jPave. Qppenheimer. fc Sivanandam I l2008l ). It also 
suppresses star formation in agreement with high-redshift 
luminosity functions ([Pave. Finlator fc Qppenheimer I 
l2006l : iMarchesini et all l2007l ) and the cosmic evol ution of 
UV luminosity at 2; ~ 4 — 7 (|Bouwens et al. IlioOTi ). Hence 
although this model is heuristic and does not describe the 
detailed physics of wind driving, it appears to regulate 
star formation in high-z galaxies in broad accord with 
observations. 

The primary simulation analysed here contains 512^ 
gas and 512^ dark matter particles, in a random cubic 
periodic volume of 96 Mpc/h (comoving) on a side with 
a gravitational softening length (i.e. spatial resolution) of 
3.75 kpc/h (comoving, Plummer e quivalent). We assurn e 
a WMAP-5 concordant cosmology (jKomatsu et al. Il2008h . 
specifically 0.^ = 0.28, J7a = 0.72, Hq = 70 km/s/Mpc, 
Q.b = 0.046, n = 0.96, and as = 0.82. This yields gas and 
dark matter particle masses of 1.2 x 10® M© and 6.1 x 10* Mq, 
respectively. The initial conditions are generated with an 
lEisenstein fc Hu I (|l999l ) power spectrum at z = 199, in the 
linear regime, and evolved to 2; = 0. We have found through 
resolution convergence tests that this is the coarsest possible 
mass resolution that adequately models the star formation 
histories of 2 ~ 2 — 4 galaxies. Given that this quarter billion 
particle simulation took several months on « 100 Harper- 
town cores, it represents the state of the art for studying rare 
large galaxies within a random cosmologically-representative 
volume. Fortunately, its volume is well-matched to current 
surveys of SMGs. 

We identify galaxies using Spline Kernel Interpolative 
DENMAX ( SKID), and halos u sing a spherical overdensity 
criterion; see lKeres et al.l (|2005l ) for more details. A galaxy's 
instantaneous star formation rate is the sum of star forma- 
tion rates of all star-forming particles in the galaxy. When we 
quote metallicities, they are star formation rate-weighted, 
since this corresponds most closely to nebular emission line 
measures of gas-phase metallicities. Resolution tests indicate 
that galaxies with > 128 star particles have well-converged 
star formation histories; for our 96 Mpc/h box, this corre- 
sponds to 7.7 X IQ^Mq. We take this as our galaxy mass 
resolution limit. Note that stella r mass functions are con- 
verged to at least half that mass (|Finlator et al. |[2006l ) , but 
we are being conservative here. In any case, this paper will be 
concerned with substantially more massive systems. There 
are 6437 galaxies at 2 = 2 in our resolved galaxy sample. 



3 SIMULATED SMGS 

The defining characteristic of SMGs is their high bolometric 
luminosity, which most likely implies a high star formation 
rate. Hence to identify SMGs in our models, we make the 
ansatz that SMGs are the most rapidly star-forming galax- 
ies in our simulated universe. In that case, we can identify 
simulated SMGs as all galaxies above a chosen SFR thresh- 
old whose number density matches the observed number 
density of SMGs. At 2 ~ 2, the observed observed num- 
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Figure 1. The stellar mass function (in number per dex per 
Mpc^) at 2 = 2 from simulations (lines) and observations of 
IMarchesini et ahl l|2009l . points). The dashed line shows results 
from our fiducial 96 Mpc//i volume, while the dotted green line 
shows a higher-resolution 48 Mpc/h volume; the overlap demon- 
strates good numerical convergence. The blue line shows the mass 
function of our simulated SMG sample. Data points are shown at 
2 < z < 3 (red) and 1.3 < z <2 (magenta). Agreement with the 
simulations is good, particularly for M > 10^^ Mq, which is key 
for this study. Note that the quoted observational errors do not 
include systematic uncertainties in the redshift determination or 
SED fitting; see iMarchesini et all l|2009l) for those. 



ber density is 1-2 x lO"'^ Mpc"^ (jChapman et al. II2OO5I : 
iTacconi et allbOOSl ). Adopting a value of 1.5 x lO"'' Mpc"^, 
the required threshold for star formation rate in our sim- 
ulations is 180 Mq/jt. None of our conclusions are signif- 
icantly altered had we chosen a SFR threshold higher or 
lower within the observationally-allowed range. Qur cho- 
sen threshold yields 41 galaxies in our simulated volume 
of 2.58 X 10^ Mpc^ We will call this our simulated SMG 
sample. 



3.1 Stellar mass function 

In Figure [1] we show the 2 = 2 stellar mass function of our 
fiduci al simulation (da s hed l ine) compared to recent data 
from IMarchesini et all (|2009l) at 2 < 2 < 3 (red points) 
and 1.3 < 2 < 2 (magenta). The agreement is reasonably 
good, particularly at the massive end that is relevant for 
this work. This is critical because the number density of 
large galaxies (which is a key barometer of SMG models) 
depends sensitively on feedback presc riptions, cosmology , 
etc. We note that t he mo dels used in iFardal et al.l (|200ir i 
and IFinlator et al. I (|2006l ) overpredicted the number densi- 
ties at the b rig ht end. We suspect that the calculation of 
iDekel et al.l |2009 ) . in which most of the gas accreted at 
the virial radius ends up forming stars, would also show 
this flaw0 Hence our current simulation, and particularly 
its momentum-driven wind feedback prescription, provides 
a more plausible platform to investigate the nature of SMGs. 

This figure also demonstrates numerical resolution con- 
vergence of our stellar masses. The green dotted line shows 



^ In addition, the simulation volume of iDekel et al. I 1I2OO9I) is 
seven times smaller than that used here, so it cannot statistically 
model a population of galaxies as rare as SMGs. 
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the stellar mass function from a higher-resolution simulation 
identical in every way to our fiducial simulation (including 
the same number of particles), except with half the box size 
(48 Mpc/h), and half the softening length (1.875 kpc/h). 
The agreement in the overlapping mass range with the 
larger- volume simulation is excellent. Overall this shows that 
our simulation produces a robust and observationally con- 
sistent set of galaxy stellar masses at z = 2. 

The simulated SMGs (blue line) occupy the most mas- 
sive end of the stellar mass distribution. They range in mass 
from lO" - 1O^^M0, and above 3 x 10^^ M© nearly every 
galaxy is identified as a SMG. Hence, SMGs are among the 
most massive galaxies at this epoch in our simulations. The 
fact that the most rapidly star forming galaxies coincide 
with the most massive galaxies is our first key result, and 
our first prediction for the nature of SMGs. 

3.2 Stellar masses 

Figure [2] (top left) shows star formation rates versus stel- 
lar mass for our simulated SMG sample (green circles) and 
for our resolved galaxy sample (red points). The simulated 
SMGs lie on the SFR— M* relation defined by the lower 
mass galaxies, i.e. these systems are not outliers as would 
be expected if their SFRs are being substantially boosted 
by mergers. In fact, there is nothing particular that distin- 
guishes SMGs from the remainder of the galaxy population, 
with the exception of their large stellar masses. 

The observational data on SMG stellar masses currently 
span a wide range, likely owing to the numerous systemat- 
ics involved with deriving stellar masses from photometric 
near-IR data. We show a selection of data from various au- 
thors in the figure. A key point to note is that the sys- 
tematic differences in stellar masses between observational 
samples are significantly larger than the range of stellar 
masses withi n any given sam ple. To illustrate this, one can 
compare the iHainline I ||2008| ) data (cyan squares) with the 
iMichalowski et all (|2009F data (blue triangles). These two 
analyses, for the most part, use exactly the same data for 
exactly the same galaxies, but employ two diffe rent algo- 
rithm s to obtain their stellar masses - yet the iHainlinel 
l|2008l ) stellar masses are on average factor of 6 lower! In 
general, other analyses have tended to find masses agreeing 
with the higher mass range; the SHA DES Lockman Hole 
galaxies analysed by iDve et al] (120081 ) are shown as the 
magenta squares, with SFRs determined fro m a subset of 
those for which 350/im data was obtained bv lCoppin et al.l 
l|2008[) to more ac c urate ly ch aracterise the du s t tem pera- 
ture. iBorvs et al. I (|2005l ) and ISwinbanlc et al.l l|2006D also 
find typical SMG masses well above IO'^^Mq. However, it 
may be that these analyses have substantially overestimated 
the stellar masses owing to contamination in the rest-frame 
H-ba nd (observed Spitzer/lRkC) by AGN (Hainlinc ct al. 
|2009| ). In contrast. [Michalowski et ai.l (|2009ll argue that the 
AGN contamination is small owing to a similar far IR-to- 
radio flux ratio in SMGs as that found in local star-forming 
galaxies. Clearly, the last word has not yet been spoken on 
the stellar masses of SMGs. 

Our simulations agree w ell with the SHADES data 
and IMichalowski et al.l l|2009l 'l analysis showing high stellar 
masses. This is a strong motivati ng factor for pref erring our 
interpretation of SMGs. If instead IHainline 1 (|2008l ') were cor- 



rect and the other determinations are biased too high, then 
our conclusions would be signiflcantly altered, and our sce- 
nario would not be as viable. Determining the stellar masses 
of SMGs to even within a factor of two would be a major 
step towards understanding their nature. 



3.3 Star formation rates 

Turning to the star formation rates, here there is clear dis- 
agreement between the simulations and the observations: the 
simulated SMGs have SFRs that are systematically too low 
by a factor of a few compared to the real SMGs. Observation- 
ally, SMGs show SFRs from ~ 400 - 2000 Mo/yr (derived 
from far-IR data), whereas the simulated SMG range from 
180 — 570 Mq/jt. Hence our star formation rates are too 
low by a factor of ~ 2 — 4. All the observational analyses 
agree for the most part, si nce they ty p ically just translate 
far-IR flux into SFR using iKennicutt I (|l998al ). though sig- 
nificant uncertainties still remain particularly owing to poor 
constraints on the dust temperature. 

At face value, the mismatch in SFRs would seem to be 
a catastrophic failure of this model; clearly these galaxies 
need some significant boost of star formation, such as might 
occur in a major merger. Our simulations do not have the 
numerical resolution to follow the detailed dynamical pro- 
cesses within merging galaxies that trigger a starburst, hence 
one might suppose that they are systematically missing such 
a strongly starbursting population. 

However, the answer is not so simple. As we pointed 
out in the introduction, hierarchical models predict that the 
major merger rate (i.e. < 3 : 1 ) of such massive ga laxies is 
around unity per Hubble time ||Guo fc White ir2008l ). Hence 
of our 41 simulated galaxies, one would expect that approx- 
imately 1 should be undergoing a major merger, assuming a 
~ 100 Myr duration for the SFR boost. Indeed, we see one 
galaxy with a significantly boosted SFR, which upon inspec- 
tion is undergoing a merger: It has a SFR of 574 Mq/jt (the 
highest among our simulated SMGs) , with a stellar mass of 
2.8 X 10^^ M0, whereas a more typical galaxy at this stel- 
lar mass has a SFR of ~ IbOA^Q /yr (we will examine this 
galaxy further in |4]). It may be that had we resolved this 
galaxy better the starburst might be stronger and be in bet- 
ter agreement with the observed SMG star formation rates. 
But in the larger picture, this cannot explain the boosted 
SFR seen in essentially all SMGs relative to the model ex- 
pectations. Hence, mergers might obtain the correct far-IR 
fiuxes but not the correct number density at these large 
masses. 

One way to reconcile all of this is if SMGs probed merg- 
ers further down the mass function. If one out of every 30 
galaxies were undergoing a merger (which is a maximal as- 
sumption, since the major merger rate falls to lower masses), 
this would imply that SMGs would have to have masses 
down to ~ 3 X 10^" M© to achieve the correct number den- 
sity. In that case, the observed large stellar masses would 
be poorly fit in the models, although this might be better 
accommodated if the stellar masses from iHainline , (,20081 ) 
are correct. This is (qualit atively) the d i lemm a faced by 
the semi-an alytic models of iBaugh et al.l (|2005l ). which as 
discussed in lSwinbank et aL I (|2008l ) match the number den- 
sity but fail to match the stellar masses. Furthermore, such 
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Figure 2. Top left: Star formation ra te versus stellar mass for the simulated S MG samp le (green circles), all other s imulated galaxies 
(red points), and observed SMGs from lHainline I l|2008l ) (cyan squares). iMichalowski et al.l |2009) (blue triangles), and lDve et al \ l|2008l l 
(magenta squares). The red line shows the median SFR— M* relation for simulated galaxies , with a Icr scatter shown by the error bars. 
The cyan line shows a fit to SFR— M* for star-forming BzK galaxies from lPaddi et all l|2007t ). with approximate la error bars (±0.3 dex) 
indicated by the dashed lines. Top right: Gas fraction /gas = Mgas/(Mgas + M,) for simulated SMGs (green) and gala xies (red), plus 
a running median (red line). Black circles show gas frac tions inferred for z ~ 2 BM/BX galaxies by lErb et al. I l|2006l ) . Blue squares 
show more direct gas fractions from CO measurements bv lTacconi et al~ 1 1I2OO8I '). Bottom left: Local galaxy density smoothed over 1 Mpc 
spheres. SMGs tend to live in denser environments, but par ticularly at lower masses there is a substantial spread. Bottom right: Oxygen 
metallicity. Data points for BM/BX galaxies are shown from lErb et al. I {2OO6). SMG metallicities tend to be higher than typical galaxies, 
following the mass-metallicity relation defined by lower-mass systems. Note that the lower-mass SMGs tend to have higher gas fractions, 
lower metallicities, and (by construction) higher star formation rates than other simulated galaxies of the same M«. 



low-mass systems would seem to require an implausibly top- 
heavy IMF to produce the SMG far-IR luminosities. 



tional data, and to gain a greater physical insight into the 
nature of this population. 



In short, our simulations broadly match the number 
densities and the stellar masses of SMGs but fail to repro- 
duce the observed SFRs by a modest factor. Other popular 
scenarios fail to match different ones of these three quanti- 
ties. It is not clear where the answer lies, but in ij5]we will 
examine how the SFR discrepancy might be symptomatic 
of a more general problem with 2; ~ 2 galaxies. For now we 
will focus on examining other properties of our simulated 
SMGs to see how our scenario fares against other observa- 



3.4 Gas fractions 

The upper right panel of Figure [5] shows the gas fractions 
(= Mgas/(Mgas + Mt)) of our simulated galaxies, with the 
subsample of SMGs indicated by the green points. Binned 
mean gas fractions from a sample of BM/BX galaxies (i.e. 
2: ~ 2 ph o tomet rically-selected galaxies) are shown from 
lErb et al~l l|2006l ). It is worth cautioning that this com- 
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parison is fairly crude; first, the simulated gas masses are 
taken to be the mass of star-forming gas in each galaxy, 
which may not correspond directly to what would be inferred 
from m olecular or Hi emission. Furthermore, the lErb et al. I 
l|2006l ) data do not come from direct measurements of gas 
masses, but rather are gas masses inferre d from the star for - 
mation rate surface density assuming the lKennicutt I (|l998l ) 
relation. Despite these caveats, the simulated galaxy popula- 
tion shows a general agreement with the trend of higher gas 
fractions in lower-mass systems. Note that this is critically 
dependent on our outflow model; other outflow models d o 
not match as well (jPave. Finlator. fc Qppenheimer Il2007l l. 

Direct measures of gas fractions at z ~ 2 — 3 are now 
becoming available thanks to i mproving millimetre i nter- 
ferometry. CO gas measures by iTacconi et al.l (j2008l ) are 
shown as blue squares in the figure. Unfortunately, once 
again the comparison is not straightforward, since the obser- 
vations include only H2-|-He, while the simulations include 
all the cold gas including Hi. Also, the observed measure- 
ments depend on the famous X factor, the conversion factor 
from CO emis s ion to molecular hydrogen mass; the data of 
iTacconi et al.l (|2008l ) is best fit by an X factor similar to 
that found in local ULIRGs, which is about 5x lower than 
that in the Milky Way disk. Still, the data generally lie in 
the range of the s imulated galaxies . The one exception is 
the "Cosmic Eye" (jSmail et al.ll2007l : ICoppin et alll2007l ) at 
Mt = 1.5 X 1O^°M0, which is the blue point weU below 
the relation at low masses; it could be a 2a outlier, or its 
M* or /gas could be misestimated, or it may have a large 
Hi reservoir. Even so, there are still simulated galaxies with 
similarly low /gas at these stellar masses. 

Although simulated SMGs broadly follow the /gas — M* 
anti-correlation seen in lower-mass systems, it is notable 
that at a given stellar mass, SMGs tend to be more gas- 
rich objects. This simply reflects the fact that SMGs are 
a SFR-selected sample, which in turn selects galaxies that 
have large gas reservoirs to fuel star formation. While the 
simulated galaxies are broadly in agreement with the data, 
more direct measures of gas mass, together with a better 
understanding of the relationship between Hi and molecu- 
lar gas, can potentially provide more robust and interesting 
constraints. 

3.5 Metallicities 

The bottom right panel of Figure [5] shows the metallic- 
ity of simulated galaxies versus M«. The SMGs lie along a 
well-defined mass-metallicity relation that broadly matches 
an extension of the relation for z ^ 2 B M/BX galax- 
ies jErb et al. II2006I : iFinlator fc Dave II2008I ). Because they 
are extremely massive, they also have quite high metallici- 
ties, typically solar or more. Note that the amplitude agree- 
ment is perhaps not so meaningful given the systematic un- 
certainties in metallicity measures of distant galaxies (e.g. 
iKewlev fc Ellisoii1l2008l ). but the relative metallicities (and 
hence the trend with M*) are likely more robust. 

At a given M, , the SMGs tend to have lower metallic- 
ities; this goes alo ng with their gas richne ss and high star 
formation rate. In IFinlator fc Dave I (|2008l ) we pointed out 
that galaxy metallicities are set by a competition between 
accretion (which dilutes the metallicity) and star formation 
(which increases it), the balance of which is modulated by 



outflows. Hence, a galaxy that has recently acquired a large 
amount of gas through an infall event will have a lower 
metallicity, as well as a higher gas content and star forma- 
tion ra±e;_riu^_preo^ in local galax- 
ies (|Ellison et al.ll2008l : iPeeples et aTl 120091 '). It wiU be in- 
teresting to see, as the data improves, whether this trend is 
also true in SMGs compared to similar mass-selected galaxy 
samples. 



3.6 Environment and clustering 

The bottom left panel of Figure [5] shows the local galaxy 
density of our simulated galaxy sample spherically averaged 
over 1 comoving Mpc, plotted versus stellar mass. The local 
density is computed using all the resolved galaxies (M* > 

7.7 X 10^ Mq). Our simulated SMGs clearly prefer to live 
in overdense environments. However, at a given stellar mass 
there is no particular trend; SMGs and non-SMGs are not 
obviously distinguished by environment. Hence the dense 
environs simply reflect the high degree of clustering expected 
for a massive galaxy population. 

Observational estimates of SMG environments are dif- 
ficult, owing to the relati vely small and n on-un iform cover- 
ages of existing surveys. ISerieant et al. I l|2008l ) found that 
SMGs tend to prefer dense environments (though not ex- 
clusively so), although t he statistics w e re onl y definitive for 
1 < z < 1.5 systems. iTamura et al.l (|2009l ) studying the 
SSA22 protocluster regions at z = 3.1 al so finds evidence 
that S MGs prefer to live in denser regions. IChapman et al. 1 
finds that SMGs appear to be more biased than 
Lyman break galaxies. The idea that SMGs tend to live 
in dense environments is qualitatively consistent with our 
model, but a quantitative comparison will have to await im- 
proved data sets, particularly over wider areas; SCUBA-2 
will help in this effort. 

Our large simulation volume allows us to make predic- 
tions for the clustering of SMGs. Figure |3] (top panel) shows 
the autocorrelation functions of simulated SMGs (solid 
black) at z = 2, along with that of our entire resolved galaxy 
sample (having a number density of 2.5 x 10~^ Mpc~^, sim- 
ilar to that of bright Lyman break galaxies) and the dark 
matter. The autocorrelation of SMGs clearly shows a larger 
amplitude than the full galaxy sample, while its slope is 
fairly similar. The clustering length of simulated SMGs is 
ro « 10/i~^Mpc, compared to 5.2ft~'^Mpc for the full re- 
solved galaxy sample. 

In the bottom panel we show the bias of SMGs relative 
to all resolved galaxies (green) and to the dark matter (blue). 
We estimate the bias by taking the square root of the ratio 
of the correlation functions. Relative to the dark matter, the 
bias of SMGs is ~ 6, which is about twice that of typical 
galaxies. In short, SMGs are highly biased and clustered 
objects in our simulation, as would be expected from their 
large stellar masses. 

Observations of SMGs support t he idea that they are 
highly clustered. iBlain et al. I l|2005h found a correlation 
length of ro ~ 7.5 ± 2.6/i"^Mpc for S MGs, suggesting that 
they are located in fairly massive halos. iFarrah et al] (|2006l ) 
determined ro ~ 6 ± l/i"^Mpc for ULIRGs at 1.5 <« < 3, 
though many of these had IR luminosities below that of 
SMGs. These determinations are somewhat lower than we 
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Figure 3. Top: Correlation function of our simulated SMGs 
(black solid), of our resolved galaxy sample (green dashed), and 
the dark matter distribution (blue dotted). Error bars on the 
SMG curve show Poisson errors in the pair counts; the Poisson 
errors for the other curves are negligibly small and are likely dom- 
inated by cosmic variance, which is not included here. Bottom: 
The bias of SMGs estimated from \/ relative to the dark 
matter (blue dotted), and of all galaxies (green dashed). SMGs 
are a highly clustered population with tq ~ 10/i~^Mpc and a bias 
of ~ 6. 



predict, although given the uncertainties, the agreement is 
not bad. 

The general prediction that simulated SMGs live in 
fairly dense and clustered environments seems to be borne 
out, at least preliminarily, in available data. But a detailed 
comparison will have to await large uniform samples of 
SMGs as is anticipated from upcoming instruments. Clus- 
tering may provide a key discriminant of SMG scenarios, 
since if SMGs are mergers of more ordinary-sized galax- 
ies, they will generally cluster as such, although a weak ef- 
fect is expected from merger bias (i.e . that m ore massive 
galaxies have a higher merger rate; e.g. IGuo fc White .,.2008i : 
IChapman et al. II2009I ). If SMGs are massive, they have little 
choice but to live in massive halos that cluster strongly. 



4 EXAMPLE SIMULATED SMGS 

We now concentrate on several individual simulated SMGs 
to study their properties in more detail. We choose four 
galaxies that span the range of SMGs from our z = 2 simu- 
lation output: 

A. The highest SFR SMG (M. = 2.8 x 10"Mq, 
SFR=574M0/yr); 

B. The most massive SMG (M* = 7.9 x lO", 



Figure 4. Top: Star formation histories in 50 Myr bins for the 
four galaxies described in text. The right edge of the histograms 
corresponds to z = 2, the adopted epoch of observation. Bottom: 
Mean metallicity of star formation as a function of time for the 
four chosen galaxies. 

SFR=334M0/yr); 

C. The median M. SMG (M. = 2.7 x 
SFR=224M0/yr); 

D. The least massive SMG (M. = 1.1 x lO", 
SFR=216A/0/yr). 

We also refer the reader to our website, 
: //luca. as . arizona. edu/~oppen/IGM/subnim.html[ 

where we show images similar to Figure [5] for all of our 
simulated SMGs. 



4.1 Star formation and enrichment histories 

Figure |4j top panel, shows the star formation histories 
(SFHs) of these four galaxies. Overall, the SFHs are rela- 
tively constant since 2 ~ 6 (t = 1 Gyr), with the larger 
stellar masses simply owing to a higher rate of quiescent SF. 
Excursions at the ~ x2 level from the time-averaged SFR 
are not uncommon. 

However, in their recent histories, i.e. in their last 
~ 200 Myr, our SMGs divide into two types, ones whose 
current SFR is comparable to their recent average (B,C) 
and ones that are experiencing a current boost in their 
SFR (A,D). This is directly related to where they lie in 
the SFR— Af, plane (Figure [21): A and D are clearly out- 
liers above the median SFR, while B and C he close to the 
median. Ifence, we can generalise these trends by looking at 
where SMGs lie relative to other galaxies in the SFR— M« 
plane. 

The careful reader will note that the final z — 2 SFRs 
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in Figure |4] are not always identical to the SFRs listed for 
these objects in ^ The difference is that the plot shows 
50 Myr-averaged star formation rates, while the numbers 
quoted above are instantaneous SFRs. For the more quies- 
cent systems, B & C, the numbers are similar but for the 
bursting systems, A & D, the instantaneous rate is non- 
trivially higher than the past 50 Myr smoothed average. 

While the scatter in SFR— Af» is fairly small (0.3 dex), 
this is still comparable to the spread in SFRs for SMGs, 
which means that a SFR-selected sample like SMGs will 
preferentially pick out galaxies like A and D that are un- 
dergoing a boost. Therefore, despite an underlying trend in 
SFR— Mt, examining only SMGs in current samples would 
likely result in a scatter plot in this plane, as shown in Fig- 
ure [21 As one probes to lower masses, SMGs more strongly 
select high-SFR outliers in SFR-M.. 

Figure |4l bottom panel, shows the growth of metallicity 
in these galaxies. All four galaxies enrich themselves to at 
least 50% solar by 2 = 6, showing that even with strong 
outflows that drive the majority of metals into the IGM, 
massive galaxies still enrich themselves quite rapidly. The 
most massive galaxy has the earliest formation epoch and 
hence the earliest metal growth, but by z ~ 4 all the galax- 
ies have similar metallicities around solar. A consequence of 
rapid enrichment is that there is also expected to be large 
amounts of dust in these galax;ies from a fairly early epoch 
(assuming that dust production tracks metallicity on a rea- 
sonably short timescale), providing the obscuration required 
to produce the large far-IR luminosities. We leave for fu- 
ture work an analysis of dust repr ocessing in these systern s 
to predict a detailed SED (as in iNaravanan et al.l l2009al ): 
this will likely require substantially higher spatial resolution 
than we can currently achieve in representative cosmological 
volumes. 

4.2 Morphologies and kinematics 

One argument in favour of a merger origin for SMGs is their 
apparently compact and disturbed morphologies. Double 
nuclei, tidal tails, and kinematics that depart strongly from 
simple spiral or elliptical g alaxies are seen in many cases (e.g. 
iTacconi et aLll2006l . [200a ). How can a model that postulates 
that these objects are oversized but otherwise normal star- 
forming galaxies be reconciled with such a scenario? The 
short answer is, at this e poch most star-fo rming galaxies are 
compact and disturbed (|Giavaliscd |2002| ) . hence such traits 
are not necessarily i ndicative of an o ngoing merger as would 
be the case locally (|Law et al.ll2007l ). 

Figure [5] shows 50 x 50 kpc (physical; about 12" at 
z = 2) maps of our four example simulated SMGs, in vari- 
ous physical quantities. Galaxy A (left panels) is undergoing 
a merger at this epoch with a smaller object about one- 
eighth its mass, seen just to the lower left. This results in 
a more centrally-concentrated distribution of gas and stars 
compared to the other systems, as well as a fairly chaotic ve- 
locity field. We note that the smaller object was originally 
significantly larger, about twice its current mass, but was 
tidally disrupted to its current size during the merger. 

The largest galaxy, B, appears the most quiescent and 
ordered of them all. It shows an extended gaseous star- 
forming centre with a large central galaxy showing ordered 
rotation. It is being fed by gas along various streams, but 



while some lumps are apparent, the majority com es in rel- 
atively smoothly. As argued in iKeres et al.l (|2009l ). that is 
typically the case for high-z galaxies. 

The other two galaxies (C,D) show more typical exam- 
ples of simulated SMGs. In both cases they show relatively 
extended gas and star-forming regions, at least compared to 
A. However, the velocity fields do not show any particular 
order, as there are a fair number of other small galaxies in 
the vicinity as well as infalling gas resulting in somewhat 
chaotic kinematics. 

The actual appearance of these galaxies in observable 
bands depends on a number of other effects. CO traces neu- 
tral gas density to some extent, modulo variations in the 
molecular gas fraction and the ratio of molecular gas to CO. 
Ongoing star formation can be traced directly in the far-IR 
continuum, although the relationship of the reradiating dust 
to the star formation may not be trivial. Radio observations 
may provide a more direct constraint, assuming that any 
contribution from a central AGN can be resolved out. The 
stellar density can be most directly traced in the rest near- 
IR, but given the amount of dust in SMGs, even those wave- 
lengths can be substantially affected by extinction and its as- 
sociated patchiness. The velocity fie ld traced in Ha with in- 
tegra l field units such as SINFONI l|F6rster Schreiber et al.l 
I2OO9I ) provides a fairly direct measure of the star forming 
gas kinematics, if not extincted. All these observations are 
possible today, but their interpretation is c omplex. Dust and 
molecular radiative transfer models as in INaravanan et aU 
(|2009bf ) are required for a proper interpretation; in a fully 
cosmological setting this is a strong challenge for numericists 
in the coming years. 

Using millimet r e inte rferometry plus other wavelength 
data. FTkcconi et al.l (|2006l ) provide a variety of evidence that 
SMGs are scaled-up versions of local major mergers seen 
as ULIRGs. Many of the arguments involve the compact- 
ness of the star formation as seen in CO lines, typically 
1 — 2 kpc, although these were high-CO transitions that typ- 
ically come from only the densest gas. On the other hand, 
there is good evidence from radio interferometry that the 
star formation in SMGs is extended over many kpc scales 
(as in our simulated SMGs in Figu re [5l) and not confined 
to the nucleus as in local ULIRGs (|Chapman et al. 1 12004 
iBiggs fc Ivisonll2008l ). Also, SMGs appear significantly less 
reddened than nuclear starburst s, which furthe r differenti- 
ates them from nearby ULIRGs l|Hainline et al.ll2009f ). This 
may argue against the major merger interpretation, since in 
mergers the star-forming gas is generally driven into a very 
compact and highly obscured nuclear starburst that in lo- 
cal ULIRGs are < 1 kpc, although simulations of SMGs 
from mergers do show larger extents of a few kpc (e.g. 
INaravanan et al]|2009bl ). Hence, current results on CO mor- 
phologies are not conclusive at to whether or not SMGs are 
major mergers or not, and it remains possible that SMGs 
could be a heterogeneous population. The kinematics are 
also not yet conclusive, as even our non-merging systems 
show disordered kinematics. A key advance in this area 
would be high-resolution imaging in the near-IR, with, e.g., 
Hubble/WFC3 or JWST, to probe the stellar distribution 
on sub-arcsec scales (modulo patchy extinction). We pre- 
dict single smooth stellar distributions in most cases with 
perhaps small companions, whereas a pre-coalescence major 
merger would clearly show two distinct large concentrations. 
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Figure 5. Simulated 2 = 2 SMGs A, B, C, D (left to right), in slices 50 kpc (physical; 6") across, showing from top to bottom the stellar 
density, gas density, star formation rate density, and line-of-sight velocity field. 



5 RELATIONSHIP TO OTHER GALAXY 
POPULATIONS 

5.1 Are the observed star formation rates 
overestimated? 

The overall properties for SMGs derived from simulations 
broadly agree with the available observations of SMGs, in- 
cluding their stellar masses, gas fractions, morphologies, and 
clustering. However, there is still one clearly discrepant is- 
sue: the star formation rates in the simulated galaxies are 
a factor of a few below that inferred from far-IR measures 



of SMG SFRs. Is this discrepancy sufficiently serious to rule 
out our scenario for SMGs, or is there a plausible explana- 
tion? 

To investigate this, we introduce the star formation ac- 
tivity parameter 

a,f = (M./SFR)/(te-lGyr) (1) 

(jPave l[2008l ). asi measures the fraction of a Hubble time 
(minus the first Gyr in which little growth happens) that 
this galaxy must form stars at its current rate to produce 
its current stellar mass. 

Figure [6] shows a^i for the simulated SMGs (green 
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Figure 6. The star formation activity parameter Ogf = 
(M,/SFR)/{tH - 1 Gyr) for the simulated SMGs (green) and 
all the simulated galaxies (red), with green sohd hues showing 
the running median for the SMGs. The cyan dashed line shows a 
Umit of SFR= 180 Mp,/y i at 2 = 2. Obse r vation s of SMGs from 
iMichalowski et al.l (120091) and iDve et al.l ||2008| ') with the SFRs 
arbitrarily lowered (and thus Q^f raised) by a factor of 3 are plot- 
ted as blue triangles and magenta squ ares, respectively , with a 
running median shown for the Michalo wski et al ] l|2009l ) sample. 
Such a rescaling places the median value of a^t in line with sim- 
ulated SMGs. 



points), which indicates that a^i ~ 0.8 for simulated galax- 
ies i n this mass range. Taking the data of iDve et al.l I2OO8I ) 
and IMichalowski et al] l|2009l ) at face value yields a typical 
Qsf « 0.2 — 0.3. If we apply an ad hoc correction factor of 
x3 to their a^f, then the results are as shown by the blue 
and magenta data points. This results in good agreement 
between the simulated and observed SMGs, with a median 
Qsf ~ 0.7 — 0.8. This shows that to reconcile our model with 
data, Af*/SFR must be overestimated in the models or un- 
derestimated in the data by about a factor of three. 

Let us first consider whether the simulations might un- 
derestimate the star formation rates by ~ x 3, perhaps owing 
to some details of the subgrid star formation prescription. 
It turns out this is not viable in a cold mode accretion- 
driven scenario, where the star formation rate is strongly 
limited by the gas supply rate, which in tur n is set pre- 
dominantly by the gravitationa l accretion rate |Keres et al.l 
I2OO5I : iDekel fc Birnboim I [2OO6I : iKeres et al.ll2009D . ThaTTs 
the reason why Oaf ~ 1 is a fairly generic result attained 
in b oth numerica l and semi-analytic galaxy formation mod- 
els (iDave I [20081 ). Hence, the model galaxies cannot form 
stars faster unless gravity brings gas in faster. Alternatively, 
we argued in ii3.3l that mergers are unlikely to be responsi- 
ble for such an underestimate, since only a small fraction of 
galaxies are undergoing a merger at any given time. 

One could equivalently match a^i by claiming that the 
observed stellar masses are underestimated by ~ x3. How- 
ever, this would then require galaxies with quite large stellar 
masses, which grossly exacerbates the problem of producing 
enough such systems by 2 ~ 2 in a hierarchical cosmology. 
Note that the recently discuss ed stellar population i ssue of 
thermally pulsating AGB stars (|Maraston et al.ll2005l ) would 
tend to lower, not raise, the inferred stellar masses, though it 



is not expected to be a large factor compared to the current 
systematic uncertainties. 

Hence for our scenario to be viable, the observed star 
formation rates in SMGs must have been overestimated by 
a (modest) factor of 3. It seems rather presumptuous to 
claim based on our models that observations are wrong, but 
a factor of 3 does not seem out of the realm of possibility, 
given the current uncertainties involved in converting vari- 
ous fluxes to star formation rates in high-2: galaxies. 

As an aside, a subtle predicted trend is that a^f is cor- 
related with M*. This is characteristic of a SFR-selected 
sample; the dashed cyan line in Figure [6] shows our adopted 
SFR limit of 180 Mg/yr. Smaller SMGs would be preferen- 
tially selected to have elevated SFRs (and hence lower a^i), 
while more massive galaxies are more quiescent and so would 
have cvsf ~ 1. This trend is also seen in the observed sam- 
ples but only weakly, probably owing to the currently large 
systematic uncertainties in determining M« and SFR. 

There are other lines of evidence that suggest that SFRs 
in 2 ~ 2 galaxi es of a l l type s have been overestimated. As 
pointed out in iDavel ([2008'), this would re concile the ob- 
served values of Qsf in BzK galaxies jPaddi et al. 2007,1 with 
the simulated Oaf. More such evidence comes from exam- 
ining the global cosmic star for mation rate density evolu- 
tion versus stellar mass growth jWilkins et al. I [20081 ): they 
showed that the time derivative of the cosmic stellar mass 
density at 2 ~ 2 is lower by ~ x3 compared to the star 
formation rate inf erred from high-mass tr acers such as Ha 
or UV (though see lReddv fc Steidel llioogi ). suggesting that 
the true SF R is lower t h an tha t obtained using conversion 
factors from iKennicutt I (|l998al ). 

Possible physical reasons for such an overestimate are 
numerous. It could simply be a calibration issue; the locally- 
calibrated conversions between far-IR, UV or radio luminosi- 
ties could be off by this factor in the case of high-2 galaxies. 
Galaxies as active as SMGs do not exist in the local Uni- 
verse, so it is difficult to test such calibrations in analogous 
systems nearby. One could imagine that the physical con- 
ditions in the ISM are systematically different, since high-2 
galaxies tend to have much higher star fo rmation rate sur- 
face densities than typical galaxies today (|Erb et al. II2OO6I : 
iTacconi et al ] |2008l ). It could also be an issue with dust tem- 
peratures, as the inferred SFR is h ighly sensitive to its as- 
sumed value (e.g. lFardal et al. 2001). Observations of 350/im 
and 850/im fluxes bv lCoppin et al.. (,200 8) indicate a typical 
dust temperature of 28 K (with substantial scatter), which 
is lower than the canonica lly-assumed 35 K from local star- 
bursts jPunne et aLllioOol ): this will be better determined 
with upcoming Herschel data. The lower temperature could 
imply as much as x2 lower SFRs. There also may be some 
far-IR light from AGN, particularly for the most luminous 
sources, although direct estimates show that it is likely to 
be sub-dominant in SMGs. Finally, SED fltting itself can be 
uncertain at a factor of ~ 2—3 level in SFR wh en only broad- 
band data are employed l|Muzzin et al. |[2009l ). although it's 
not clear that this would provide a systematic shift. It is not 
a stretch to think that one or more of these factors could 
lead to a factor of 3 overestimate in the SFRs. 

Another perhaps related possibility is that the IMF 
could be different in high-2 galaxies. This has been offered 
as a solution to various nagging quandaries in galaxy and 
stellar evolution, such as the colour and luminosity evolu- 
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tion of early-type galaxies (|van Dokkum I [2OO8I ) . the abun- 
dance of_caxbon enhanced metal-poor stars in the Milky 
Way (jTumlinson Il2p07[). and the evolution of the cosrnic ste l- 
lar mass assembly (jPardal et al.ll2007l : IWilkins et al. II2OO8I ). 
One so lution i nvoke d to solve all these problems, as pointed 
out in iDave I (|2008l ). would be that the IMF is somewhat 
more top-heavy or bottom-light at higher redshifts, such 
that SF Rs would be o verestimated by ~ x3 at 2 ~ 2. 
Recently. iMeurer et al.l (, 2009 ) suggested that nearby galax- 
ies with higher SFR surface densities have more top-heavy 
IMFs, which adds some local credibility (albeit controver- 
sial) to the idea of the IMF perhaps being more top-heavy 
(or bottom-light) at early epochs. 

The level of IMF var iation required in o ur scenario is 
far milder than that in the lBaugh et al.l l|2005l ) semi-analytic 
models of SMGs, and it is not easy to constrain directly even 
in local galaxies let alone at z ~ 2. It is, for instance, quite 
consistent with dire ct constraints on the IMF of SMGs by 
iTacconi et al.l (|2008l V IMF variations would also have other 
implications, such as more metals (and hence dust) being 
produced for a given stellar mass, which could have impor- 
tant effiects_on_ch]stt and inferred SFRs (see 
e.g. iBaugh et al.ll2005h : a thorough analysis of all its effects 
is beyond the scope of this work. It must be pointed out 
that there is no clear evidence for IMF variations in any 
context, so this possibility should be considered with some 
skepticism. 

For our scenario, it is irrelevant whether the cause of the 
purported SFR overestimates is a different IMF or some cal- 
ibration effect. The main point is that if the SFRs have been 
systematically overestimated by a factor of 3, then this re- 
moves the main clearly discrepant aspect of our model when 
compared to current observations of SMGs. Such an overes- 
timate may not be unique to SMGs, but may be reflective 
of a systematic difference in all z ~ 2 star-forming galaxies 
relative to the local ones. 



5.2 Duration of the SMG phase 

SMGs in our scenario are forming stars relatively quies- 
cently, as seen in Figure [l] Hence one expects that SMGs 
remain identifiable as SMGs for a long time. This contrasts 
with the short duty cycle of SMGs in a merger-based sce- 
nario. Inferred duty cycles for SMGs tend to be short, but 
these are mostly based on gas consumption timescales. In 
our scenario, such timescales are not relevant because the 
IGM supplies gas at a rate comparable to the star forma- 
tion rate. Therefore, our SMGs don't "use up" their gas 
but merely consume it at the rate at which it is being sup- 
plied. The key point is that simulations generically predict 
remarkably large gas supply rates in high- 2: galaxies. 

To get a handle on the duration of the SMG phase, we 
identify the progenitors and descendants oi z = 2 SMGs 
in our simulations at z = 2.5 and z = 1.5, i.e. 0.67 Gyr 
prior and 1 Gyr later. In Figure [7] we show the location in 
SFR— A/, of 2: = 2 SMGs (green points) at these earlier and 
later epochs. For reference, the horizontal line demarcates 
the SFR threshold needed to maintain a number density of 
1.5 X 10~^ Mpc~'^; the values are slightly different than at 
2 = 2, namely 130 Mq/jt at 2 = 2.5 and 190 M^/yr at 
2 = 1.5. 

At both 2 = 2.5 and 2 — 1.5, slightly more than half 
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Figure 7. SFR vs. M« at 2: = 2.5 (top panel) and z = 1.5 (bottom 
panel). Descendants/progenitors of z = 2 simulated SMGs are 
shown as green points. Blue squares denote SMGs that would 
be identified based on a SFR threshold at those redshifts; the 
SFR thresholds are indicated by the dotted lines. More than half 
the z = 2 SMGs are identified with similarly-selected SMGs at 
z = 2.5 and z = 1.5. 



the 2 = 2 SMGs are among the galaxies that would be iden- 
tified as SMGs based on an analogous SFR threshold. The 
remainder fall just below the SFR threshold. In general, the 
2 = 2 SMGs at these other epochs are fairly randomly dis- 
tributed in SFR amongst the general galaxy population at a 
given Mt . This is in line with our scenario that SMGs have 
a fairly constant SFR, with minor (~ x 2) excursions owing 
to stochastic accretion events. 

Expressed in terms of a duty cycle, one might say that 
the SMGs in our simulations have a duty cycle of ~ 50%, 
since at any given epoch the high-SFR criterion picks out 
about half the massive galaxy population. This duty cycle 
depends on mass; the lower-mass SMGs tend to be greater 
outliers, and hence have smaller duty cycles. Of course, the 
concept of a duty cycle may not be exactly appropriate for 
these systems, given that they do not become dramatically 
less active when they are not in the SMG phase, but rather 
simply fall slightly below what would be identified in an 
SFR-selected sample. 

One implication of this relatively high duty cycle is that 
the number density and clustering of SMGs should inde- 
pendently translate into roughly the same halo mass. This 
would indicate that a large fraction of halos at the high- 
mass end are occupied by SMGs. This is in contrast to a 
scenario where SMGs are major mergers, in which the small 
duty cycle of SMGs in a merger phase would imply a low 
fraction of massive halos containing an SMG. With wider 
and more uniform samples, this prediction can be tested 
relatively straightforwardly, which we expect will strongly 
discriminate between the different SMG scenarios. 
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Figure 8. Histogram of halo masses of simulated SMGs at 2 = 2 
(red solid line), and the halos in which their descendants live at 
z = 1 (green dotted) and 2 = (blue dashed). The median values 
at each epoch are indicated at the top of the figure. 



5.3 SMG descendants today 

Where do the descendants of SMGs live today? In our sim- 
ulation we can follow our simulated SMG sample to 2: = 0, 
where we can directly identify the halos in which they reside. 

Figure [8] shows a histogram of the halo masses of simu- 
lated SMGs at 2 = 2 (red solid line) , and the halos in which 
their descendants live at 2 = 1 (green dotted) and 2 = 
(blue dashed). At each epoch, these are among the most 
massive halos in the simulation. From 2 = 2 — >■ 0, the typi- 
cal halo mass grows from 1.7 x lO^M© to 1.1 x 10^* Mq. The 
latter is a group-sized halo. Hence, we can infer that SMGs 
are typically the brightest group galaxies in the process of 
active formation at 2 ~ 2. This is also consistent with their 
number density and large clustering length. 

Another issue is that SMGs cannot continue to form 
stars at their 2 ~ 2 rates all the way to 2 = 0, since no 
galaxies approaching those bolometric luminosities are seen 
locally, and their stellar masses would be far too large. Our 
simulated SMGs match the 2 ~ 2 observed stellar mass func- 
tion (Figure [T]). How much longer can they sustain these 
rates and still be consistent with lower-2 mass functions? 

To answer this, let us make the ansatz that SMGs 
turn into the most massive galaxies today. Clearly this is 
not strictly true, but given that the simulated SMGs are 
among the most massive galaxies at 2 ~ 2 living in the 
most massive halos, it may not be a bad approximation. 
In any case it provides a upper limit to how much SMGs 
can grow. Examining the stellar mass functions versus red- 
shift of|Marchosinij5t j|L (2009.), it is clear that there is not 
much room for growth at the massive end. For SMG num- 
ber densities of 10~^ Mpc~^, the stellar mass has grown by 
only 0.2 dex since 2 ~ 2. Hence even if one were to ignore 
dry merging, this only allows for a growth of 60% in stel- 
lar mass in the last 10 Gyr. If the typical stellar mass is 
M, = 1O"-^M0, then this is a growth of ^ 1.2 x IO^Mq . 
Taking a median SFR of « 4OOM0/yr (|Coppin et al.ll2008h . 
this means a typical SMG can only continue forming stars at 
this rate for at most 300 Myr, and perhaps much less given 



that massive galaxies also grow by dry merging. Therefore, 
one would be observing SMGs at a fairly special time, when 
the star formation has to shut off within <^ O.lin. If, how- 
ever, the SFRs were lowered by a factor of three, the SMGs 
could conceivably continue to form stars at their observed 
rates for the better part of a Gyr before quenching, and 
hence SMGs would not be being observed at a particularly 
special time. The latter interpretation is more consistent 
with our scenario that SMGs are not typically undergoing a 
spectacularly life-changing event. We note that there would 
still be some phenomenon required to quench star formation 
in massive galaxies to produce the red and dead systems seen 
today, but perhaps it does not have to be directly associated 
with the SMG phase. 



6 CONCLUSIONS 

We investigate the physical properties of the most rapidly 
star-forming galaxies at 2 = 2, constrained to match the 
observed number density of sub-millimetre galaxies, in a 
ACDM hydrodynamic simulation of galaxy formation. We 
examine their stellar masses, star formation rates, metal- 
licities, environments, clustering, star formation histories, 
morphologies, and kinematics. We find that the simulated 
high-SFR galaxies are a good match to the observed SMGs 
in terms of their stellar masses, number densities, and clus- 
tering, which is a first for a cosmologically-situated galaxy 
formation model. 

However, there is a significant discrepancy between the 
predicted and observed star formation rates: The simulated 
SMG's SFRs are lower by a factor of ~ 3 compared to that 
inferred for SMGs from their observed far-IR fiuxes. We ar- 
gue that it is plausible that the SFRs of SMGs have been 
overestimated by this factor. There is growing evidence that 
such a modest overestimate would reconcile various obser- 
vations related to star formation and stellar mass growth 
at 2 ~ 2. Given that high- 2 galaxies are in many ways sys- 
tematically different than present-day ones (e.g. they have 
much higher star formation rate surface densities), it is pos- 
sible that the conversion factors calibrated to local galaxies 
may be incorrect at high-2. One way to alter such calibra- 
tions would be to invoke an IMF with more massive stars 
relative to low-mass ones. This IMF variation is within cur- 
rently allowed constraints, and is far milder than the ex - 
tremely top-heavy IMF proposed in the lBaugh et al.l (|2005| ) 
semi-analytic models. However, there are many other possi- 
ble reasons why the SFR calibration could be off for SMGs. 

We note that the agreement with the observed stellar 
masses depends with which observational sample one com- 
pares. Most analyses seem to favour large stellar masses for 
the SMGs, typically > 10^^ M©, but some obtain signifi- 
cantly lower stellar masses, even using identical data. This 
illustrates the current uncertainties in determining stellar 
masses from photometric data, especially in the presence of 
large amounts of extinction and possible AGN contamina- 
tion, in a regime that is poorly calibrate d by local data. I f the 
lower stellar masses as inferred by e.g. iHainlinel (|2008l ) are 
confirmed, our model would be significantly less preferred, as 
it would require an order of magnitude overestimate of the 
observed star formation rates to reconcile our simulations 
with such data. Furthermore, it would alleviate one of the 
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difficulties witti tlie major merger model, namely that large 
major mergers are too rare, since these lower-mass systems 
are more common. However, it is still true that the merger 
model would have to explain the enormous far-lR fluxes, 
which seem to require quite massive gas-rich progenitors or 
else a dramatic change in the IMF. 

If our scenario is correct, the most significant impli- 
cation would be that SMGs are typically not large major 
mergers of gas-rich disks, as is canonically believed. Hence, 
they are in this sense not analogs of local ULIRGs, in spite of 
their high bolometric fluxes. They are instead, to first order, 
super-sized versions of ordinary star forming galaxies, sitting 
at the high-mass end of the star formation rate-stellar mass 
relation. They are being supplied with gas from minor merg- 
ers and smooth cold accretion at a rate comparable to their 
SFR. The SMG phase is, therefore, not a particularly dra- 
matic event in the life of a massive galaxy, although their 
star formation must eventually be quenched by some mecha- 
nism for these systems to become brightest group ellipticals 
today. 

It is worth noting that our simulations do not pro- 
duce a population of passively evolving galaxies as observed, 
likely because we do not explicitly include any physical pro- 
cess that fully quenches star fo rmation, such as AGN feed- 
back fe.g. ldi Matteo et al.ll2008|). Passively evolving galaxies 
are seen at z ^ 2 (e.g. iKriek et al. I l2008h. and at the mas- 
sive end perhaps as many as half the galaxies are passively 
evolvi ng, or at least forming stars at a reduced rate for their 
mass (jPapovich et al. II2OO6I ). If we declared half our mas- 
sive galaxies to be passively evolving, then to reproduce the 
number density of SMGs, we would have to go to a some- 
what lower SFR threshold of 120 Mq/jt (vs. 180 Mq/jt). 
Hence, to match the observed SFRs we would then require 
a reduction by a factor of ~ 4 in SFR instead by a factor of 
~ 3 owing to some systematic, which is still not implausi- 
ble. Of course a full galaxy formation model should account 
for passively evolving systems, but the fact that our current 
simulations do not is unlikely to affect our basic conclusions 
here. 

Although SMGs are not particularly distinguished from 
the more typical galaxy population other than by their larger 
stellar masses, they are still effectively a star formation rate- 
selected sample. Therefore, we predict that there will be mild 
trends to higher SFRs, lower metallicities, and higher gas 
fractions than for non-SMGs at comparable stellar masses. 
In other words, the specific star formation rate of SMGs will 
be higher than the average at their stellar mass, and this will 
be particularly true for lower-mass SMGs owing to the scat- 
ter in the SFR— Af* relation. Stated in terms of the star for- 
mation activity parameter Osf = (Af«/SFR)/(fH — 1 Gyr), 
one expects that SMGs at higher masses will have a higher 
asf , since they will be less likely to have a temporarily ele- 
vated SFR. If this trend is confirmed in future studies (which 
will require much-improved measures of SFR and Af*) it 
would provide further support for our overall scenario. 

Our model predicts further trends that should become 
evident with larger and deeper SMG samples: (i) Lower-Af» 
SMGs will appear preferentially more disturbed or merger- 
like, while high-A^« systems will appear more quiescent; (ii) 
Going a factor of a few deeper in far-IR luminosity should 
start to suggest a trend in SFR— A/, for SMGs; (iii) This 
trend should smoothly join onto the trend of SFR— Af, 



seen in optically-selected galaxies; (iv) SMGs should be 
highly clustered, and the clustering strength should increase 
with SMG luminosity; (v) The occupancy fraction of SMGs 
within massive halos, as might be quantified by compar- 
ing the clustering and number densities of SMGs, should be 
fairly high. These trends are likely subtle, so testing them 
will require improved modeling along with improved data. 
We particularly emphasize the importance of accurate stellar 
mass and clustering measures to provide optimal discrimi- 
nation between models. Naively, one might expect that the 
above trends would not exist or would be weaker in a ma- 
jor merger-dominated scenario for SMGs, but this must be 
assessed within a full hierarchical model. 

All the current scenarios for SMGs, i.e. the large ma- 
jor merger scenario, the modest merger scenario with very 
top-heavy IMF, and our super-sized star-forming galaxy sce- 
nario, each have their own difficulties. Our main point is that 
the latter scenario is at least as consistent with the avail- 
able observations of SMGs as all the others. Hence, there 
should be some doubt as to whether or not SMGs are sys- 
tems caught in the act of a major merger. Even so, major 
mergers of large galaxies almost certainly occur, and based 
upon the simple number density estimates we presented ear- 
lier (SJII, it might be expected that perhaps one or two in 
ten SMGs are indeed such systems. In fact, this may be the 
only way to produce the most extremely luminous SMGs, 
e.g. with Lib. > 10^"^ Lq. Therefore, we speculate that per- 
haps SMGs are a heterogeneous population, and that the 
very most extreme systems are indeed large major merg- 
ers, while the more common modestly-luminous SMGs are 
super-sized star formers. 

SMGs are a key phase in the evolution of massive galax- 
ies, possibly the most active phase of formation for present- 
day passive spheroids. Assembling a fully self-consistent 
cosmological framework for their formation and evolution 
will require much effort on both theoretical and observa- 
tional fronts. Fortunately, there will be substantial obser- 
vational advances forthcoming in long-wavelength capabili- 
ties to trace star formation, dust continuum emission, and 
molecular lines of SMGs out to high-z. These will be com- 
plemented by better multiwavelength data to constrain their 
stellar masses, AGN content, clustering, and metallicities. 
Together with rapidly improving models of SMGs within a 
cosmological framework, there is growing hope that the na- 
ture of these enigmatic systems and their importance in our 
overall picture of galaxy formation will soon be revealed. 
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